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PREFACE 

The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing is a multiyear program of research, development, evaluation, and appli- 
cation of aerospace remote sensing for agricultural resources, which began in 
fiscal year 1980. This program is a cooperative effort of the U.S. Department 
of Agriculture, the National Aeronautics and Space Administration, the National 
Oceanic and Atmospheric Administration (U.S. Department of Commerce), the 
Agency for International Development (U.S. Department of State), and the 
U.S. Department of the Interior. 

The work which is the subject of this document was performed by the Earth 
Resources Applications Division, Space and Life Sciences Directorate, 

Lyndon B. Johnson Space Center, National Aeronautics and Space Administration 
and Lockheed Engineering and Management Services Company, Inc. The tasks 
performed by Lockheed Engineering and Management Services Company, Inc., were 
accomplished under Contract NAS 9-15800. Funding for this study was provided 
by the National Oceanic and Atmospheric Administration under the Foreign 
Commodity Production Forecasting (FCPF) project of the AgRISTARS program. 
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1, INTRODUCTION 


The Agriculture and Resources Inventory Surveys Through Aerospace Remote Sensing 
(AgRISTARS) program is a multiyear program of research, development, evaluation, 
and application of aerospace remote sensing. The goal of the program is to 
determine the usefulness, cost, and extent to which remotely sensed data can be 
integrated into existing or projected U.S. Department of Agriculture (USDA) 
systems to improve the objectivity, reliability, timeliness, and adequacy of 
information required to carry out USDA missions. 

Major efforts in the AgRISTARS Foreign Commodity Production Forecasting (FCPF) 
project and its historical predecessors have been directed toward the estima- 
tion of crop acreages from spectral data obtained by the Landsat multispectral 
scanner. Over time, a number of cultural and environmental parameters have 
been identified which cause difficulty in the interpretive process. One such 
parameter is the crop moisture condition at the time the spectral data are 
collected. 

A number of approaches to the monitoring of crop moisture conditions are being 
developed and evaluated in the context of the AgRISTARS program. Unfortunately, 
a common feature of these relatively sophisticated technologies is the require- 
ment for data which are currently available only in a research mode. The FCPF 
project is particularly constrained by the need to use widely available data 
that are at least potentially available in real time. In this paper, an interim 
solution to the problem of estimation of crop moisture conditions is presented. 
This solution is the form of an indicator of relative agricultural moisture 
conditions and is designed for international application within the realistic 
expectations for availability of meteorological data obtained from group 
observations. 

The index is based upon a moisture balance between precipitation (gains) and 
evaporation (loss) for a uniform medium. A key element of this index is the 
estimation of evaporation from maximum and minimum air temperatures. 
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The following section contains a brief review of the technical literature on 
evaporation and evapotranspiratlon estimation. Details of the development of 
the evaporation model are given In section 3. The moisture Index Is presented 
In section 4. 


2. BACKGROUND AND REVIEW OF THE LITERATURE 

Evapotranspiratlon is crucial In the life cycle of a plant. Its primary func- 
tion is to maintain a favorable daylight temperature In the plant, that Is, to 
keep the leaves a few degrees cooler than the environment under adequate mois- 
ture conditions. The consequences of low moisture are complex: photosynthesis 
decreases due to the higher leaf temperatures caused by the decrease In 
transpiration as well as the lower supply of raw materials (water and carbon 
dioxide). 

The crop modeler finds that evapotranspiratlon measurements are even scarcer 
than those of evaporation. They may take sone consolation from Chang (ref. 1) 
who states that "a relationship — good enough for many agricultural purposes 
usually exists between evaporation and evapotranspiratlon for a period of a day 
or longer." He goes on to cite numerous examples of the superiority of Class A 
pan evaporation data over other methods of estimating evapotranspiratlon. 

However, one is still left with the lack of Class A pan data. Since both 
evaporation and evapotranspiratlon are controlled by the atmospheric demand for 
moisture (air temperature, vapor pressure deficit, windspeed, and net 
radiation), many estimation techniques use combinations of these elements. 

In 1948, Pennwn (ref. 2) derived a semi -empirical equation to estimate evapora- 
tion from a free water surface. His technique required daily mean air tempera- 
ture, dew point temperature, mean windspeed, and net radiation. The air and 
dew point temperatures were transformed to their corresponding saturation vapor 
pressures to approximate the atmospheric demand. 
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The modified Jensen-Haise method (ref. 3) estimates potential evapotranspiration 
front daily mean temperature, net radiation, and station-specific constants. One 
of the constants, the air temperature coefficient, is derived from the satura- 
tion vapor pressures associated with the normal maximum and minimum temperatures 
during the month with the greatest mean air temperature. 

Simpler methods have also been devised which require only air temperature meas- 
urements. One such evapotranspiration method is that of Thornthwaite (ref. 4) 
which uses mean monthly temperatures and an empirical annual heat index. More 
recently, several techniques to estimate pan evaporation from air temperatures 
alone have been developed. Griffiths (ref. 5) and his colleagues Moe (ref. 6) 
and Miller (ref. 7) tested and expanded a technique of estimating pan 
evaporation rates from air temperature. Their equation was of the form 

EP » ag + aJX (1) 


where 

EP » the amount of monthly total evaporation from the pan in inches 
TX * the monthly mean of daily maximum temperatures in degrees Fahrenheit 
3(5 and a^ are linear regression coefficients with Oq in inches and aj in 
inches per degrees Fahrenheit 

Their approach determined the two coefficients for each station in a network. 
Next, the coefficients were plotted on maps which were analyzed for isolines of 
these coefficients. Such maps could be used to estimate the coefficients for a 
hypothetical pan at any station observing TX. Moe (ref. 6) created such maps 
for Texas, and Miller (ref. 7) expanded the area to include most of the 
southern United States. 

The patterns of the isolines on these maps have a reasonable physical and areal 
interpretation. The theoretical basis for the formula, the special character- 
istics of the naps, and the sources of errors in the model were first inter- 
pretea and analyzed by Trenchard (ref. 8). He found that much of the error in 
the naps for Texas was due to the seasonal variability of vdndspeeds and that 
this could be reduced with a simple correction term based upon the month of the 
year. Again, maps of the coefficients were generated for interpolation. 
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The iflodels previously listed rely on the high correlation of pan evaporation 
with air temperature. This correlation varies with locati-^ primarily because 
of the geographic dependence of the relationships between emperature and the 
active factors governing evaporation. In addition, while a station's temp- 
erature may represent a region the size of a crop reporting district, the sta- 
tion's evaporation measurements may not be !ue to advection effects. However, 
evaporation models based on solely air temperature are still somewhat station- 
specific and not entirely satisfactory for large area crop modeling. It is 
believed that the statior dependency of such models may be decreased by using 
the saturation vapor pressures associated with maximum and minimum temperatures 
instead of the temperatures themselves. The use of vapor pressure equivalents 
of the temperatures helps to account for the nonlinearity of the vapor pressure 
function in the evaporation process, the variability of atmospheric moisture, 
and a portion of the influence of advection and ventilation. Thus, they are 
included in the model which is derived in the following section. 

3. MODEL DEVELOPMENT 


3.1 DATA 

The data for the development and test of this model consist of 1104 monthly 
mean maximum and minimum temperatures (degrees Fahrenheit), windspeeds (miles 
per hour), and monthly total evaporation (inches) collected at 26 weather sta 
tions throughout the Great Plains. The locations of these stations are given 
in table 1 and figure 1. The stations are part of the cooperative network 
maintained by the Environmental Data Information Services of the National 
Oceanic and Atmospheric Administration (EDIS/NOAA). 


The saturation vapor pressure in millibars corresponding to each temperature 
was calculated with ths following formula. 


Vapor « 6.11 X exp 


-176204.2621 + 5597.607915 x T - 2.850772636 x T^ 
l'2F4Tff.2+ ■273"xl 


where T is the temperature in degrees Fahrenheit. 
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148259 Kan. Tuttle Creek Dam 39°15' 96 “36* 1057 4 .72 

148648 Kan. Webster Dam 39°25' 99 “25* 1863 2 .36 

322188 N. Dak. Dickinson 46“53' 102“48’ 2585 22 3.99 

Experimental Station 


















Figure 1.- Map showing station locations in the U.S. Great Plains. 
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Half (552) of the observations for each year at each station were placed in the 
development data set. The remaining observations were reserved for testing, 

3.2 PROCEDURE 

The evaporation model was derived through regression analysis. Much of the 
statistical computations were performed with procedures available in the 
Statistical Analysis System (SAS) computer system (ref. 9) installed on the 
Laboratory for Applications of Remote Sensing (LARS) IBM 370 computer. 

The first step in model development was to i.tudy the correlation table of the 
weather variables (table 2). Pan e-- operation is most highly correlated with 
the vapor-pressure corresponding to maximum temperature, then with maximum 
temperature, minimum temperature, and minimum vapor pressure. It is least 
correlated with wind. 

As expected, the temperatures are extremely well correlated with their corres- 
ponding vapor pressures [correlation coefficient (r) = 0.98]. This suggests 
fhat, in a linear regression model, only one variable from each temperature- 
vapor pressure pair needs to be included. For example, maximum temperature and 
maximum vapor pressure provide similar information to a model. To include both 
of these variables would be redundant. 

The correlations between maximum and minimum air temperatures and maximum and 
minimum water temperatures are also great (r > 0.90). These four variables are 
similar, while the wind variable differs from them. In modeling terms, just 
one of the first four should be sufficient before adding wind. 

To obtain the model, we used forward selection of variables in the stepwise 
regression procedure available on the SAS. Through this process, the model is 
built one variable at a time. The next variable which enters the model offers 
a significant F-statistic. After each iteration, the F-statistics for all 
variables in the model are checked for significance. A variable is deleted 
from the model if its F-statistic is not significant. The procedure continues 
until either all variables have been included or the next variable selected for 
entry is deleted at the previous iteration. 
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TABLE 2.- CORRELATION MATRIX FOR THE VflSIABLES IN THE DATA SET 
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When all variables were submitted to stepwise regression, the best model included 
the following three variables; maximum vapor pressure, wind, and minimum vapor 
pressure. (Statistics and parameter estimates are given in figure 2.) Because 
wind speed data are less widely available than temperature data, the wind vari- 
able is excluded and the stepwise fitting is recalculated. This time all four of 
the remaining variables were included in the following order; maximum vapor 
pressure, minimum vapor pressure, maximum temperature, and minimum temperature 
(figure 3), It is interesting to note that the last two variables imijrove the 
model only slightly, possibly due to their high correlation with the first two 
variables. 

The comparison of the results of the stepwise regression suggests that minimum 
vapor pressure should be included in a model without the wind variable, although 
it does not improve the second model as much as the addition of the wind vari- 
able improves the first model. The coefficient of determination (R^) improved, 
and the F-va1ue remained significant at the 0.001 level. Also, the root mean 
squared error (RMSE) decreased and the bias improved. For these reasons and 
theoretical reasons to be discussed in the next section, it was believed that 
minimum vapor pressure should be retained as a model variable. 

3.3 RESULTS 

The final model estimates monthly total evaporation EP in inches from vapor 
pressures corresponding to the mean monthly maximum and minimum temperatures 
(TX and TN, respectively). 

EP ^ 0.2163 + 0.3473 Vapor(TX) - 0.2644 Vapor(TN) (3) 

where Vapor is the vapor pressure function given in equation (2). With the 
development set, equation (3) predicts pan evaporation with a RMSE of 1.72, a 
bias of 0.0012, and a correlation coefficient of 0.8232. The model performs 
equally well with the test data set and given as the RMSE of 1.71, shows a bias 
of 0.08 and an r of 0.8228. 

The use of the saturation vapor pressure function over water as a transforma- 
tion of the air temperature values is merely an effort to make the model 
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CIPI • 

23,00762169 


OF 

SUM OF squares 

MEAN SQUARE 

F 

PROU>F 

SS81”™ 

TOTAL 

59§ 
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5092.92969568 


S77.39 

0.0001 


H value 

STO ERROR 
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F 

PHOHXF 

INTFRCFPT 

VAPOR.MX 

VAPORiMN 

!:5U1K|! 

-n,269A9009 



’Sl:n 

l:mi 

5TFP 3 VARIAftlP 

MIN ENTERED 

R SQUARE » 

0.67995966 

C(M) ■ 

21.12586297 


OF 

SUM OF SQUARES 

mean SQUARE 

F 

PROtlJf 

TOTAL 

5a3 

551 

mmmi 


388.09 

O.OOOl 


n VALUE 

STO error 

TYPE II SS 

F 

PHOO>F 

INTFRCEPT 

VAPOR.MX 

VAPOR^MN 

MIN 


0.01793177 
0.08322399 
0. 03731511 

'"fillHii 


0.0001 


STEP * 

PggPgSSION 

TOT*L 


VARIAPU FNTCHEn 
OF 


intfrcfpt 

vapor^mx 

VAPOR.MN 

MAX 

MIN 


A 

•iAT 

SSI 


H VAUUE 

«.noni<J?A« 

ntAAA2AA3H 

-o.ATinArSR 

■8:?45{2?|s 


R square ■ 

SUM OF SQUARES 
vTl) ERROR 


0.69021980 C(P> 

MEAN SQUARF 


TYPE 


03 


SS 


S.OOOOOOOO 

F PR0«>F 

309.69 O.OnOl 


Slil 


F 

8S 


PROlUF 


n. ooo 

o. nno 
0.000 
0.000 


NO OTHfR VARIABIFS MFT IMF 0.5000 SIGNIFICANCE LfVEL FOR ENTRY INTO THE MODEL. 


Code: 

1. Thu source of vaHatlon regression. 

2 . Thu source of variation error. 

3. The source of variation total, 

4. Degrees of freudom. 

5. Sums of sguares. 

6. Mean squares. 

7. f-valuc, Mhich Is the ratio of the regression wan square to the error wan square. 

8. The significance probability of the F-value. 

9. R*, the square of the wUlpIo correlation coefficient. 

10, The names of the Independent variables Included In tbe model. 

11, The corresponding estimated regression coefficients. 

12, The TyiH! II sum of squares for each variable. 

13, F-vaiues and significance probabilities associated with the Type II sums of squares. 


Figure 3.- Stepwise regression procedure for dependent variable Pan 

without the wind variable. 
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variables similar to those used in more sophist i,!;ated models requiring both 
saturation vapor pressure values for air and vapor pressure values for the 
water surface. An example of this use is given by Penman (ref. 2). The 
function itself is nonlinear over the normal temperature range of agricultural 
activity, a property which should not be Ignored in the evaporation process. 

The inclusion of minimum vapor pressure is based upon the empirical relationship 
of minimum temperature with dew point. With reduced advection in the hours 
before dawn, minimum air temperature is usually limited by the dew point. If we 
assume this to be true, minimum temperature becomes significant in the applica- 
tion of the evaporation model to both moist and dry climatic regimes, where 
maximum temperature alone may not be a distinguishing factor. The negative sign 
for the regression coefficient of this variable in the model is probably more 
than coincidence and agrees well with other models which use the vapor pressure 
at the dew point, e.g., the modified Jensen-Haise method (ref. 3). 

The omission of any ventilation or advection term from the model is a signifi- 
cant source of error. Wind confounds the relationship of minimum temperature 
and dew point, usually producing a positive difference between the two and a 
serious overestimate of atmospheric moisture. On the other hand, advective 
effects are frequently reflected in higher daily maximum temperatures (the 
positive component of the model) so that some compensation is possible under 
windy conditions. 


4. SPONGE, A GENERALIZED MOISTURE INDICATOR 

In part, the evaporation algorithm was developed in order to assess crop mois- 
ture status. In a general sense, crop moisture is the balance of precipitation 
and evapotranspiration over tin>e. 

The crop moisture index (CMI) (refs. 10 and 11) Is currently used In many 
AgRISTARS projects to assess moisture conditions. It includes a E-layer soil 
water model and potential evapotranspiration calculated with Thornthwaite's 
method. Some of its requirements restrict its use to regions for which 
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long-term average data as well as current precipitation and temperature data 
are available. An additional restriction Inherent In the CMI Is that It Is an 
indicator of regional moisture for periods of at least a week and may give 
unreliable results when applied to a single station on a dally basis. 

The evaporation algorithm of equation (3) allowed us to develop a simple 
moisture Indicator with a sound physical basis that used common meteorological 
variables, was suitable over a broad range of climates, and was applicable to a 
single station. The result was named "sponge." Figure 4 shows a conceptual 
Illustration of sponge. 

Sponge is described as a simple medium with 8 Inches of water-holding capacity 
which Is initialized half-full of water on January 1*. Each day, in accordance 
with the hydrologic cycle, water Is added to the medium from precipitation and 
lost through evaporation. Precipitation (both liquid and frozen) is added at 
the full amount until the layer Is saturated. It is this sponge-like behavior 
which gives the variable its name. Any additional precipitation is assumed to 
be lost as run-off or drainage. Evaporation occurs at a fraction of the Class 
A pan rate, the exact proportion being the ratio of the current contents to the 
total capacity of the sponge. Either actual or estimated evaporation pan 
values may be used. The daily contents of the sponge are defined as 

S. » S^.i + - (E^ X S^.i/CAP) (4) 

where 

- sponge contents on day i, in inches 

P^ a precipitation on day i, in inches 

Ej a actual or estimated pan evaporation in inches on day i 
CAP * sponge capacity in inches 

and 0 < < CAP. 


*Aliernati vely, the final value of the previous year may be used as an initial 
value, and the capacity may be varied for a particular region. 
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When evaporation pan nieasureiiients are not available, they may be estimated with 
equation (3), with a divisor of 30 days to convert the evaporation function to 
a daily value. 

Si - - (EP(TX^, TN^) X S^.i/CAP x 30) (5) 


where 

EP * pan evaporation function (equation 3) 

TX^ » maximum ten^e nature on day 1 
TN^- ■- minimum temperature on day i 

Because of its simple data requirements (daily precipitation and evaporation 
estimated from maximum and minimum tentperatures), the sponge can be calculated 
at any temperature-precipitation observation station. 

5. SUMMARY 

An attempt has been made to align meteorological observations with crop condi- 
tions by developing a method of estimating Class A pan evaporation (a surrogate 
for evapotranspi ration) from air temperatures. A regression model of monthly 
evaporation totals from 26 Class A pans in the U.S. Great Plains accounted for 
two-thirds of the variance of this variable. Ventilation (wind) accounted for 
most of the remaining variance but was not included so that the model could 
have more general application. 

The first suggested use for these evaporation estimates has been to define a 
hypothetical medium (sponge) and to propose a simple budget of precipitation 
and evaporation as an indicator of meteorological moisture stress. The simple 
form and minimal data requirements for this new variable make it an ideal 
candidate for investigation and application to various AgRISTARS projects which 
use meteorological data. 
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